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Lattice QCD determination of states with spin | or higher in the spectrum of nucleons 
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Energies for excited isospin i states that include the nucleon are computed using quenched, 
anisotropic lattices. Baryon interpolating field operators that are used include nonlocal operators 
that provide G2 irreducible representations of the octahedral group. The decomposition of spin | 
or higher states is realized for the first time in a lattice QCD calculation. We observe patterns of 
degenerate energies in the irreducible representations of the octahedral group that correspond to 
the subduction of the continuum spin | or higher. 



PACS numbers: 12.38.Gc 21.10.Dr 

The theoretical determination of the spectrum of 
baryon resonances from the fundamental quark and gluon 
degrees of freedom is an important goal for lattice QCD. 
To date there have been studies of experimental ground 
state energies for different baryons [1-3] but only a few 
results for excited state energies have been reported. [4-9] 
No clear determination of states with spin | or higher has 
been reported because nonlocal operators have not been 
used. In this work we find degenerate energies that oc- 
cur in irreducible representations of the octahedral group 
corresponding to the subduction of the continuum spin | 
or higher. This provides the first lattice QCD calculation 
that realizes the decomposition of spins greater than |. 

Lattice correlation functions do not correspond to def- 
inite values of total angular momentum. However, they 
do correspond to definite irreducible representations (ir- 
reps) of the octahedral group when the source and sink 
operators transform accordingly. There are six double- 
valued irreps of the octahedral group: three for even- 
parity that are labeled with a g subscript (gerade) and 
three for odd-parity that are labeled with a u subscript 
(ungerade). They are: Gig, Hg,G2g,Giu, Hu and G2u- 

Continuum values of total angular momenta are real- 
ized in lattice simulations by patterns of degenerate en- 
ergies in the continuum limit that match the patterns for 
the subduction of spin J to the double-valued irreps of 
the octahedral group. These patterns are shown in Ta- 
ble L For example, a state in one of the G2 irreps is a 
signal for the subduction of continuum spin | or higher. 



TABLE L The number of occurrences of double-valued irrep 
A of the octahedral group for different values of continuum 
spin J. 
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For spin | , there must be partner states in the H and G2 
irreps that would be degenerate in the continuum limit. 
For spin ^, there must be partner states in the Gi, H 
and G2 irreps and for spin | there must be one part- 
ner in the Gi irrep and two in the H irrep. The partner 
states should be degenerate in the continuum limit where 
lattice spacing a ^ 0. 

This paper reports on work to determine the pattern 
of low-lying states in the isospin ^ spectrum. We carry 
out an analysis in quenched lattice QCD using a mod- 
erately large number of three-quark operators. Work is 
in progress to use a very large number of operators. [12] 
Smeared quark and gluon fields are used. [10, 11] Smear- 
ing reduces the couplings to short wavelength fluctua- 
tions of the theory and provides cleaner determinations 
of effective energies. This is important when a large array 
of interpolating field operators is used in order to imple- 
ment the correlation-matrix method of Refs. [13, 14]. 

In order to determine suitable operators, including the 
nonlocal ones that are required for the G2 irreps, we have 
developed sets of baryon operators that transform ac- 
cording to irreducible representations using an analyti- 
cal method based on appropriate Clebsch-Gordan coeffi- 
cients for the octahedral group. [15] An alternative, auto- 
mated procedure that is developed in Ref. [16] provides 
very large sets of operators. Both methods provide equiv- 
alent results. In this work, wc use for positive parity the 
three-quark operators defined in Tables VI, VII, and X 
of ReL [15]. These comprise a complete set of quasilo- 
cal operators plus the simplest set of nonlocal operators 
that have one quark displaced relative to the other two 
quarks. 

Negative-parity operators are obtained by applying the 
charge-conjugation transformation to the positive-parity 
operators. A three-quark operator that transforms ac- 
cording to irrep A and row A of the octahedral group is 
related by charge conjugation to an operator that trans- 
forms according to irrep Ac and row Ac and has opposite 
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parity, i.e., 

Cfii^'^Ct = Csf ^)ct = -B^'^^*. (1) 

If A is Gig,G2g or Hg, then Ac is Giu,G2u or re- 
spectively, and vice versa. The row labels are related by 
Xc = + ^ — ^, where is the dimension of the irrep: 
dci =2,dH= 4, = 2. 

A Hcrmitian matrix of correlation functions for irrep A 
and row A is obtained by including a 74 matrix for each 

quark field in the source operator b\!^^\o, 0). This gives 
Cit^\t) = Pi^^ Y.i^\T{Br\.M''\o,0))\0), 

X 

(2) 

where the factor pj^^^ is ±1 and it arises from evaluat- 
ing the action of the 74 matrices in the source operator. 
Subscripts k and k' refer to different embeddings of irrep 
A. The relation between matrices of correlation functions 
with different parities is [16], 

cit^it) = -vtcitr'^^*iT-t). (3) 

where % = 1 for periodic boundary conditions and 
rjt = —1 for anti-periodic boundary conditions. Apply- 
ing charge conjugation and time-reversal to a correlation 
function produces one for the opposite parity because of 
PCT symmetry. We use a sufficiently long time extent, 
T, so that this symmetry provides two independent cor- 
relation functions for each parity. 

Analyses of excited state energies are based upon the 
matrices of correlation fimctions of Eq. (2) . In order to 
extract the spectrum of energies from the matrix of corre- 
lation functions, we first average over rows because they 
provide equivalent results owing to octahedral symme- 
try. We then solve the following generalized eigenvalue 
equation, 

E ^it^ (^)-i' ^ *o) = to) E ^It^ ito)v^^ *o), 

k' k' 

(4) 

where superscript n labels the eigenstates. The symbol 
(jW indicates that the appropriate average over rows 
has been performed. The reference time tg in Eq. (4) is 
taken near the source time t = in order to have sig- 
nificant contributions from excited states. The principal 
eigenvalues a^^^ {t, to) are related to the energy En by [17] 

a(") {t, to) e-^"(*-*°) (1 + 0(e-l^^l*)) , (5) 

where SE is the difference between £"„ and the next clos- 
est energy. We have determined effective energies En by 
fitting the generalized eigenvalues to the leading term of 
Eq. (5). 

Methods that improve the signal-to-noise ratio are very 
important for extracting excited state energies from lat- 
tice QCD simulations. The correlation-matrix method 
optimizes the operators so as to achieve early plateaus 



of the effective energies and the use of PCT symmetry 
increases the statistics. In addition, we use anisotropic 
lattices with temporal lattice spacing at one-third of the 
spatial lattice spacing a^. The finer spacing increases the 
number of time-slices for extraction of energies before the 
signal for a high-energy state decays to the level of the 
noise. 

In this work anisotropic lattices with two different vol- 
umes are used: 239 gauge field configurations are used 
for a 16'^ X 64 lattice and 167 configurations are used 
for a 24^ x 64 lattice. Because of the use of PCT sym- 
metry, the size of the statistical ensemble is effectively 
doubled. Gauge-field configurations are generated using 
the anisotropic, unimproved Wilson gauge action [18-21] 
in the quenched approximation with /3 = 6.1 and the 
pion mass is 490 MeV. For both lattices, the temporal 
lattice spacing corresponds to a^^ = 6.05(1) Gev [8] as 
determined from the string tension. Calculations are per- 
formed using the Chroma software. [22] Further details 
of the action and lattice parameters will be presented 
elsewhere. 

We have extracted energies for isospin ^ and | chan- 
nels by diagonalizing matrices of corrc;lation functions 
formed from three-quark operators that share the same 
octahedral symmetry. Starting with a large number of 
operators, we first eliminate operators that have less in- 
fiuence on the determination of the effective energies and 
the eigenvectors v^"^ of the low- lying states. Once a set 
of good operators is obtained, we form final matrices of 
correlation functions, diagonalize them and extract ef- 
fective energies. This procedure yields solid results for 
energies of low-lying states and is more efficient than di- 
agonalizing matrices of the largest dimension. We have 
obtained 17 energies for nucleonic states and 11 energies 
for delta baryon states. In this paper, we focus on the 
nucleon channel and the determination of spin. 

Our lattice results for the spectrum are shown in Fig. 1. 
For positive parity, the nucleon state is very well deter- 
mined at lattice energy atE = 0.193(3). It is an isolated 
Gig state corresponding to spin i. Well above the nu- 
cleon is a group of three degenerate states (within errors) 
in the Gig, Hg and G2g irreps. The results on both lat- 
tice volumes are essentially the same, although the lowest 
G2g state in the smaller volume has a larger error. 

Energies of G2 states are the most difficult to extract 
becaiisc the plateaus are relatively short-lived, signal-to- 
noise ratios are large, and the number of operators that 
we have used for the matrices of correlation functions is 
limited. The limitation on the number of operators is be- 
cause we restrict the operators used in this work to ones 
that can be constructed from one-link displacements. In 
work that is in progress, more varied types of operators 
are used to obtain very larger numbers of G2 operators, 
which is essential for determining states with higher en- 
ergies. 

Figure 2 shows effective masses obtained for the lowest 

energy Gig, Hg and G2g states. The horizontal lines indi- 
cate the time ranges used to fit the principal eigenvalues 
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FIG. 1: The energies obtained for each symmetry channel of 
isospin I baryons are shown based on the 24^ x 64 lattice 
data. The scale on the left side shows energies in lattice units 
and the scale on the right side shows energies in GeV. Errors 
are indicated by the vertical size of the box. 
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FIG. 2: Effective energies for the lowest states in Gig, Hg and 
G29 irreps. Time ranges used for and energies obtained from 
fits to the principal eigenvalues are shown by horizontal lines. 



Q!'^"-*(t, io) by a single exponential form. The resulting 
mean energy value and the error range are indicated by 
the lines. 

Although there are substantial discretization errors 
with the quark action that is used, and they could con- 
tribute differently in the different irreps, clear patterns in 
the degeneracies emerge. Focusing on the group of three 
positive-parity states near lattice energy atE — 0.36 in 
Fig. 1, two interpretations are possible. A.) The group 
consists of a spin i state and a spin | state that are 
accidentally degenerate. In this case the G\g state cor- 
responds to spin i and the Hg and G2g partner states 



correspond to the subduction of spin |. B.) The group 
consists of a single state with the degenerate G\g, Hg 
and G2g partner states corresponding to the subduction 
of spin |. 

In the physical spectrum of positive-parity nucleon res- 
onances, the lowest excited state, iV(1440, ^~'^), lies be- 
low all negative parity states. We do not find a signal for 
a positive-parity excitation that has lower energy than 
the negative-parity excitations at this quark mass. The 
next two excited nucleon states are essentially degener- 
ate, namely, iV(1680, f"*") and iV(1710, i^), each with a 
width of about 100 MeV. Spin | states occur only at sig- 
nificantly higher energy (1990 MeV). Primarily because 
of the absence of spin ^ in the low-lying spectrum, in- 
terpretation A.) of our lattice results is more consistent 
with the physical pattern of energies and spins. In the 
absence of G2 operators, previous lattice studies have as- 
sumed that states obtained with H irrep operators cor- 
respond to continuum spin |. [4, 6, 7, 9]. With the G2 
operators, we find low-lying states in irrep H that are 
consistent with the subduction of spin | or higher. 

In the negative-parity spectrum, we also obtain essen- 
tially the same results for both lattice volumes. The 
three lowest states shown on the right half of Fig. 1 
are unambiguously identified as follows: the lowest Giu 
state corresponds to spin ^ and the lowest two Hu states 
correspond to distinct spin | states. Above these is a 
group of three states with roughly the same lattice en- 
ergy: atE 0.33 (within errors). Again there are two 
possible interpretations. C.) The group consists of a spin 
i state in Giu that is accidentally degenerate with a spin 
I state, the latter having degenerate partner states in iJ„ 
and G2u- D.) The group consists a spin | state having 
degenerate partner states in Gi„, Hu and G2u- 

The pattern of low-lying physical states starts with 
7V(|",1520), iV(i",1535) and 7V(f",1700). These 
should show up as distinct Hu, Giu and Hu states on 
the lattice, in agreement with the three lowest negative- 
parity states in Fig. 1. The next physical states include 
7V(i~, 1680) and 7V(|~, 1675), which essentially are de- 
generate. They should show up as degenerate Giu, Hu 
and G2U states on the lattice. This pattern of spins is 
consistent with interpretation C.) of the lattice states at 
lattice energy atE 0.33. The pattern of energies of the 
physical states has ,1675) a little lower in energy 

than A^(| , 1700), but the lattice results at lattice spac- 
ing O.IF place the spin | state above the spin | state. 
Study of the continuum limit of the lattice spectrum is 
required in order to resolve these issues. 

Because the minimum spin that is contained in the 
G2 irrep is |, we have found strong evidence for spin 
I or higher for both parities in our lattice spectra for 
isospin i. We also have found evidence for degenerate 
partner states corresponding to the subduction of spin 
I or higher to the octahedral irreps. The lattice results 
for the low-lying excited states of isospin i provide the 
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correct number of octahedral states for the subduction 
of the spins of the low-lying physical states. These re- 
sults are significant but they are based on the quenched 
approximation and a 490 MeV pion mass. Similar calcu- 
lations in full QCD with several pion masses and several 
lattice spacings are planned and work is under way to 
calculate the anisotropic gauge configurations. 
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